Nuclear and nucleon transitions of the H di-baryon 



m 
o 
o 

(N 

< 



> : 
r» . 
m ■ 

oo : 
o . 

m ■ 

o ; 

^ : 

a: 
i 

D 

^ ■ 



X 



Glennys R. Farrar and Gabrijela Zaharijas 
Center for Cosmology and Particle Physics 
New York University, NY, NY 10003, USA 

We consider 3 types of processes pertinent to the phenomenology of an H di-baryon: conversion 
of two A's in a doubly-strange hypernucleus to an H, decay of the H to two baryons, and - if 
the H is light enough - conversion of two nucleons in a nucleus to an H. We compute the spatial 
wavefunction overlap using the Isgur-Karl and Bethe-Goldstone wavefunctions, and treat the weak 
interactions phenomenologically. The observation of A decays from doubly-strange hypernuclei puts 
a constraint on the H wavefunction which is plausibly satisfied. In this case the H is very long-lived 
as we calculate. An absolutely stable H is not excluded at present. SuperK can provide valuable 
limits. 



I. INTRODUCTION 

The most symmetric color-spin representation of six 
quarks (uuddss) is called the H dibaryon. It is flavor 
singlet with strangeness -2, charge 0, and spin-isospin- 
parity I(J P ) = 0(0 + ). In 1977 Jaffe calculated its mass 
[1] to be about 2150 MeV in the MIT-bag model and thus 
predicted it would be a strong-interaction-stable bound 
state, since decay to two A particles would not be kine- 
matically allowed. Since then its mass has been esti- 
mated in many different models, with results lying in the 
range 1 - 2.3 GeV. On the experimental side, there have 
been many inconclusive or unsuccessful attempts to pro- 
duce and detect it. See [2] for a review. 

The purpose of this paper is to study several processes 
involving the H which are phenomenologically important 
if it exists: conversion of two A's in a doubly-strange hy- 
pernucleus to an H, decay of the H to two baryons, and - 
if the H is light enough - conversion of two nucleons in a 
nucleus to an H. The amplitudes for these processes de- 
pend on the spatial wavefunction overlap of two baryons 
and an H. We are particularly interested in the possibil- 
ity that the H is tightly bound and that it has a mass less 
than m^v + m\. In the case, as we shall see, its lifetime 
is longer than the age of the Universe. 

If the H is tightly bound, it would be expected to be 
spatially compact. Hadron sizes vary considerably, for a 
number of reasons. For instance the nucleon is more than 
twice as large as the pion, with charge radius rjy = 0.88 
fm compared to rv = 0.38 fm. Lattice and instanton- 
liquid studies can account for these diverse radii and fur- 
ther predict that the glueball is even more tightly bound: 
r<3 w (1/4 — 1/6) r/v [3]. If the analogy suggested in ref. 
[4] between H, A1405 and glueball is correct, it would 
suggest m ~ ''g^1/4 r/v- The above size relationships 
make sense: the nucleon's large size is due to the low 
mass of the pion which forms an extended cloud around 
it, while the H and glueball do not couple to pions, due 
to parity and flavor conservation, are thus are small com- 
pared to the nucleon. In the absence of an unqucnchcd, 
high-resolution lattice QCD calculation capable of a re- 
liable determination of the H mass and size, we will con- 
sider all values of m# and take rn/r^ = 1// as a pa- 



rameter, with / in the range 2-6. For a more detailed 
discussion of the motivation and properties of a stable or 
long-lived H and a review of experimental constraints on 
such an H, see ref. [5] . 

In this paper we calculate the lifetime for decay of the 
H to various final states, and we consider two types of ex- 
perimental constraints on the transition of two baryons 
to an H in a nucleus, Abb — * A' H X. To estimate the 
rates for these processes requires calculating the overlap 
of initial and final quark wavefunctions. We model that 
overlap using an Isgur-Karl harmonic oscillator model for 
the baryons and H, and the Bethe-Goldstone wavefunc- 
tion for a nucleus. The results depend on r^/rn and the 
nuclear hard core radius. 

Experiments observing single A decays from dou- 
ble A hypernuclei Aaa[6, 7] indicate that t(Aa\ — > 
A' H X) > 10~ 10 sec. Our calculations show that adequate 
suppression of T{Aaa —> A' H X) requires r#<l/2 r^r, 
consistent with expectations, Thus an H with mass m# < 
2mA can still be viable in spite of the observation of 
double-A hypernuclei, as also found in ref. [8]. 

We calculate the lifetime of the H, in three qualita- 
tively distinct mass ranges, under the assumption that 
the conditions to satisfy the constraints from double-A 
hypernuclei are met. The ranges are m# < mjv + JTiA> 
in which H decay is a doubly- weak AS = 2 process, 
tojv + tua < ma < 2mA, in which the H can decay by a 
normal weak interaction, and tuh > 2mA, in which the 
H is strong-interaction unstable. The H lifetime, in these 
ranges respectively, is a fewx 10 12 yr, about a month, and 
- 10~ 9 sec. 

Finally, if run ^ 2m^f , nuclei are unstable to AS* = —2 
weak decays converting two nucleons to an H. In this case, 
the stability of nuclei is a more stringent constraint than 
the double-A hypernuclear observations, but our results 
show that nuclear stability bounds can also be satisfied if 
the H is sufficiently compact: Th ^ 1/4 risr depending on 
mass and nuclear hard core radius, although this option 
is vulnerable to experimental exclusion by SuperK. 

This paper is organized as follows. In section II wc 
describe in greater detail the two types of experimental 
constraints on the conversion of baryons to an H in a 
nucleus. In section III we setup the theoretical appara- 
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tus to calculate the wavefunction overlap between H and 
two baryons. We determine the weak interaction matrix 
elements phenomenologically in section IV. Lifetimes for 
various processes are computed in sections VB and VI. 
The results are reviewed and conclusions are summarized 
in section VII. 



II. EXPERIMENTAL CONSTRAINTS 
A. Double A hyper-nucleus detection 

There are five experiments which have reported posi- 
tive results in the search for single A decays from double 
A hypernuclei. We will describe them briefly. The three 
early emulsion based experiments [9-11] suffer from am- 
biguities in the particle identification, and therefore we 
do not consider them further. In the latest emulsion ex- 
periment at KEK [7], an event has been observed which 
is interpreted with good confidence as the sequential de- 
cay of He AA emitted from a S~ hyperon nuclear capture 
at rest. The binding energy of the double A system is ob- 
tained in this experiment to be Baa = I -01 ± 0.2 MeV, 
in significant disagreement with the results of previous 
emulsion experiments, finding -Baa ~ 4.5 MeV. 

The BNL experiment [6] used the (K~,K + ) reaction 
on a Be 9 target to produce S=-2 nuclei. That experi- 
ment detected pion pairs coming from the same vertex 
in the Be target. Each pion in a pair indicates one unit 
of strangeness change from the (presumably) di-A sys- 
tem. Observed peaks in the two pion spectrum have 
been interpreted as corresponding to two kinds of de- 
cay events. The pion kinetic energies in those peaks are 
(114,133) MeV and (104,114) MeV. The first peak can 
be understood as two independent single A decays from 
AA nuclei. The energies of the second peak do not corre- 
spond to known single A decay energies in hyper-nuclei 
of interest. The proposed explanation [6] is that they are 
pions from the decay of the double A system, through 
some specific He resonance. The required resonance has 
not yet been observed experimentally, but its existence 
is considered plausible. This experiment does not suf- 
fer from low statistics or inherent ambiguities, and one 
of the measured peaks in the two pion spectrum sug- 
gests observation of consecutive weak decays of a double 
A hyper- nucleus. The binding energy of the double A 
system -Baa could not be determined in this experiment. 

The KEK and BNL experiments are generally accepted 
to demonstrate quite conclusively, in two different tech- 
niques, the observation of A decays from double A hyper- 
nuclei. Therefore ta aa ->a' h x cannot be much less than 
w 10~ 10 s. (To give a more precise limit on ta aa ^a' h x 
requires a detailed analysis by the experimental teams, 
taking into account the number of hypernuclei produced, 
the number of observed A decays, the acceptance, and 
so on.) As will be seen below, this constraint is readily 
satisfied if the H is compact: rg < 1/2 rjv- 



B. Stability of nuclei 

There are a number of possible reactions by which two 
nucleons can convert to an H in a nucleus if that is kinc- 
matically allowed (rag <2m;v). The initial nucleons are 
most likely to be pn or nn in a relative s-wave, because 
in other cases the Coulomb barrier or relative orbital an- 
gular momentum suppresses the overlap of the nucleons 
at short distances which is necessary to produce the H. 
If run <2mAr — nm^ 1 , the final state can be Hir + or 
Hit and n— 1 pions with total charge 0. For nin > 1740 
MeV, the most important reactions are pn — > He + u e or 
the radiative-doubly-weak reaction nn — ► H"f. 

The best experiments to place a limit on the stability of 
nuclei are proton decay experiments. Super Kamiokandc 
(SuperK), can place the most stringent constraint due to 
its large mass. It is a water Cerenkov detector with a 
22.5 kiloton fiducial mass, corresponding to 8 10 32 oxy- 
gen nuclei. SuperK is sensitive to proton decay events in 
over 40 proton decay channels [12]. Since the signatures 
for the transition of two nucleons to the H are substan- 
tially different from the monitored transitions, a specific 
analysis by SuperK is needed to place a limit. We will 
discuss the order-of-magnitude of the limits which can be 
anticipated. 

Detection is easiest if the H is light enough to be pro- 
duced with a 7r + or 7r°. The efficiency of SuperK to 
detect neutral pions, in the energy range of interest (KE 
~ — 300 MeV), is around 70 percent. In the case that 
a 7r + is emitted, it can charge exchange to a 7r° within 
the detector, or be directly detected as a non-showering 
muon-like particle with similar efficiency. More difficult 
is the most interesting mass range tuh > 1740 MeV, for 
which the dominant channel pn — > He + v gives an elec- 
tron with E ~ (2m N - m H )/2<70 MeV. The chan- 
nel nn — ► H-f, whose rate is smaller by a factor of or- 
der a, would give a monochromatic photon with energy 
(2m N - m H ) < 100 MeV. 

We can estimate SuperK 's probable sensitivity as 
follows. The ultimate background comes primar- 
ily from atmospheric neutrino interactions, vN — ► 
N'(e,n), vN -> N'{e,n) + nn and vN —> uN' + nn, 
for which the event rate is about 100 per kton-yr. With- 
out a strikingly distinct signature, it would be difficult to 
detect a signal rate significantly smaller than this, which 
would imply SuperK might be able to achieve a sensitiv- 
ity of order T ijVN ^ x > fewlO 29 yr. Since the H pro- 
duction signature is not more favorable than the signa- 
tures for proton decay, the SuperK limit on ta nn —>a' h x 
can at best be 0.1r p , where 0.1 is the ratio of Oxygen 
nuclei to protons in water. Detailed study of the spec- 
trum of the background is needed to make a more precise 



Throughout, we use this shorthand for the more precise inequal- 
ity run < vtia ~ m A' ~ m x where mx is the minimum invariant 
mass of the final decay products. 



3 



statement. We can get a lower limit on the SuperK life- 
time limit by noting that the SuperK trigger rate is a few 
Hz[12], putting an immediate limit to^h+x > fcwlO 25 
yr, assuming the decays trigger SuperK. 

SuperK limits will apply to specific decay channels, but 
other experiments potentially establish limits on the rate 
at which nucleons in a nucleus convert to an H which 
are independent of the H production reaction. These 
experiments place weaker constraints on this rate due 
to their smaller size, but they are of interest because 
in principle they measure the stability of nuclei directly. 
Among those cited in ref. [13], only the experiment by 
Flerov et. al.[14] could in principle be sensitive to tran- 
sitions of two nucleons to the H. It searched for decay 
products from Th 232 , above the Th natural decay mode 
background of 4.7 MeV a particles, emitted at the rate 
T Q = 0.7 10 _10 yr _1 . Cuts to remove the severe back- 
ground of 4.7 MeV en's may or may not remove events 
with production of an H. Unfortunately ref. [14] does 
not discuss these cuts or the experimental sensitivity in 
detail. An attempt to correspond with the experimental 
group, to determine whether their results are applicable 
to the H, was unsuccessful. If applicable, it would estab- 
lish that the lifetime T Th 232^ H+x > 10 21 yr. 

Better channel independent limits on N and NN de- 
cays in nuclei have been established recently, as summa- 
rized in ref. [15]. Among them, searches for the radioac- 
tive decay of isotopes created as a result of NN decays 
of a parent nucleus yield the most stringent constraints. 
This method was first exploited in the DAMA liquid Xe 
detector [16]. BOREXINO has recently improved these 
results [15] using their prototype detector, the Counting 
Test Facility (CTF) with parent nuclei C 12 , C 13 and O 16 . 
The signal in these experiments is the beta and gamma 
radiation in a specified energy range associated with de- 
excitation of a daughter nucleus created by decay of 
outer-shell nucleons in the parent nucleus. They obtain 
the limits t pp > 5 10 25 yr and r nn > 4.9 1 25 yr. How- 
ever H production requires overlap of the nucleon wave- 
functions at short distances and is therefore suppressed 
for outer shell nucleons, severely reducing the utility of 
these limits. Since the SuperK limits will probably be 
much better, we do not attempt to estimate the degree 
of suppression at this time. 

Another approach could be useful if for some reason the 
direct SuperK search is foiled. Ref. [17] places a limit 
on the lifetime of a bound neutron, r„ > 4.9 1 26 yr, by 
searching for 7's with energy E 7 = 19 — 50 MeV in the 
Kamiokande detector. The idea is that after the decay 
of a neutron in oxygen the de-excitation of O 15 proceeds 
by emission of 7's in the given energy range. The back- 
ground is especially low for 7's of these energies, since 
atmospheric neutrino events produce 7's above 100 MeV. 
In our case, some of the photons in the de-excitation pro- 
cess after conversion of pn to H, would be expected to fall 
in this energy window. 



III. OVERLAP OF H AND TWO BARYONS 

We wish to calculate the amplitudes for a variety of 
processes, some of which require one or more weak in- 
teractions to change strange quarks into light quarks. 
By working in pole approximation, we factor the prob- 
lem into an H-baryon-baryon wavefunction overlap times 
a weak interaction matrix clement between strange and 
non-strange baryons, which will be estimated in the next 
section. For instance, the matrix element for the transi- 
tion of two nucleons in a nucleus A to an H and nucleus 
A', Ann — * A' H X, is calculated in the A A pole approxi- 
mation, as the product of matrix elements for two subpro- 
cesses: a transition matrix element for formation of the H 
from the AA system in the nucleus, lA^ljAAj^ff Xi times 
the amplitude for a weak doubly-strangeness-changing 
transition, \M\nn->aa- We ignore mass differences be- 
tween light and strange quarks and thus the spatial wave- 
functions of all octet baryons are the same. In this section 
we are concerned with the dynamics of the process and 
we suppress spin-flavor indices. 



A. Isgur-Karl Model and generalization to the H 

The Isgur-Karl (IK) non-relativistic harmonic oscilla- 
tor quark model[18-20] was designed to reproduce the 
masses of the observed resonances and it has proved to 
be successful in calculating baryon decay rates [19]. In 
the IK model, the quarks in a baryon are described by 
the Hamiltonian 



where we have neglected constituent quark mass differ- 
ences. The wave function of baryons can then be written 
in terms of the relative positions of quarks and the cen- 
ter of mass motion is factored out. The relative wave 
function in this model is [19, 20] 



*s(ri,r2,r 3 ) = N B exp 



4^ 



(2) 



where Nb is the normalization factor, as 



v / < r %> 



and < r% > is the baryon mean charge radius 
squared. Changing variables to 

.. fi - T2 r n + f 2 - 2r"5 



V2 



A 



V6 



(3) 



reduces the wave function to two independent harmonic 
oscillators. In the ground state 



*s(p,A) 



exp 



V + A 2 ) 



(4) 



One of the deficiencies of the IK model is that the 
value of the a b parameter needed to reproduce the mass 
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splittings of lowest lying 5 and | baryons, as = 
0.406 fm - , corresponds to a mean charge radius squared 
for the proton of < r\ h >= = 0.49 fm. This is dis- 
tinctly smaller than the experimental value of 0.86 fm. 
Our results depend on the choice oi 0.3 and therefore we 
also report results using as = 0.221 fm" 1 which repro- 
duces the observed charge radius at the expense of the 
mass-splittings. 

Another concern is the applicability of the non- 
relativistic IK model in describing quark systems, espe- 
cially in the case of the tightly bound H. With th /tn = 
1//, the quark momenta in the H are « / times higher 
than in the nucleon, which makes the non-relativistic ap- 
proach more questionable than in the case of nucleons. 
Nevertheless we adopt the IK model because it offers a 
tractable way of obtaining a quantitative estimate of the 
effect of the small size of the H on the transition rate, 
and there is no other alternative available at this time. 

We fix the wave function for the H particle starting 
from the same Hamiltonian (1), but generalized to a six 
quark system. For the relative motion part this gives 



VP h = Nh exp 



i<3 



(5) 



The space part of the matrix element of interest, 
< A' H \A A \ >, is given by the integral 

6 

J 11^^(1,2,3)^(4,5,6)^(1,2,3,4,5,6). (6) 



Therefore it is useful to choose variables for the H wave- 
function as follows, replacing 



B. Bruecker-Bethe-Goldstone Nuclear 
Wavefunction 

To describe two A's or nucleons in a nucleus we use 
solutions of the Bruecker-Bethe-Goldstone equation de- 
scribing the interaction of a pair of fermions in an in- 
dependent pair approximation; see, e.g., [21]. The solu- 
tion of the Schrodinger equation for two fermions in the 
Fermi sea interacting through a potential v{x\, x 2 ) takes 
the form 



*(1,2) = J= C 



iPRcM 



V(a) 



(10) 



where Rcm and a are defined as in (8). The first factor 
contains the center-of-mass motion and the second is the 
internal wave function of the interacting pair. tp(a) is 
a solution of the Bethe Goldstone equation (cq (36.15) 
in [21]) which is simply the Schrodinger equation for two 
interacting fermions in a Fermi gas, where the Pauli prin- 
ciple forbids the appearance of intermediate states that 
are already occupied by other fermions. Both wave func- 
tions are normalized so that the space integral of the wave 
function modulus squared equals one. In the application 
of this equation to nuclear matter, the interaction of each 
particle from the pair with all particles in nuclei through 
an effective single particle potential is included, in the in- 
dependent pair approximation known as Brueckcr theory 
(see eq (41.1) and (41.5) in [21]). 

We are interested in s-wave solutions to the above 
equation since they are the ones that penetrate to small 
relative distances. Following [21], an s-wave solution of 
the internal wave function is sought in the form 



V(a) 



u(a) 



(11) 



n, r 2 , r 3 , r 4 , r 5 , r 6 -> p a , A a , p , A , a, R CM (7) 

where ff 1 ^ and A a ( b ) are defined as in eq (3), with a(b) 
referring to coordinates 1, 2, 3 (4, 5, 6). (Since we are ig- 
noring the flavor-spin part of the wavefunction, we can 
consider the six quarks as distinguishable and not worry 
about fermi statistics at this stage.) We also define the 
center-of-mass position and the separation, a, between 
initial baryons a and b: 



Rc 



M 



Rcm + Rcm -> _ Ba Bb 

- , a — n CM n CM . 



(8) 



Using these variables, the H ground state wave function 
becomes 



/o \ 3/4 / v. 15/2 



(9) 



r oei, , ^2 -2 -2 -.2 3 
x cxp[ — ^(p a +X a +p b +X b +2 3 )]- 

As for the 3-quark system, ajj 



which simplifies the Bethe Goldstone equation to 
d 2 



(~dx 2 ~ +k ) u ( a ) = v ( a ) u ( a )- J X(a,y)v(y)u(y)dy (12) 

where v(a) is the single particle potential in the effective- 
mass approximation, and the kernel x( a >2/) is given by 



x(a,y) = - 

7T 



sin fcp(a — y) smkp(a + y) 
a-y a+y 



(13) 



For the interaction potential between two nucleons in 
a nucleus we choose a hard core potential for the fol- 
lowing reasons. The two particle potential in a nucleus 
is poorly known at short distances. Measurements (the 
observed deuteron form factors, the sums of longitudinal 
response of light nuclei,...) only constrain two- nucleon 
potentials and the wave functions they predict at inter- 
nucleon distances larger than 0.7 fm [22]. The Bcthe- 
Goldstone equation can be solved analytically when a 
hard-core potential is used. While the hard-core form is 
surely only approximate, it is useful for our purposes be- 
cause it enables us to isolate the sensitivity of the results 
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to the short-distance behavior of the wavcfunction. We 
stress again, that more "realistic" wavefunctions are in 
fact not experimentally constrained for distances below 
0.7 fm. Rather, their form at short distance is chosen for 
technical convenience or aesthetics. 

Using the hard core potential, the s-wave BG wave- 
function is 



N u(k F a) 
ly BG k F a 





for 
for 



a > 
a < 



k F 
c 

k F 



(14) 



where -S- is the hard core radius 

k F 



Expressions for u can 
be found in [21], eq. (41.31). The normalization fac- 
tor Nbg is fixed setting the integral of \iPbg\ 2 over the 
volume of the nucleus equal to one. The function u van- 
ishes at the hard core surface by construction and then 
rapidly approaches the unperturbed value, crossing over 
that value at the so called "healing distance". At large 
relative distances and when the size of the normalization 
volume is large compared to the hard core radius, u(a)/a 
approaches a plane wave and the normalization factor 
Nbg reduces to the value 1 / W, as 



i>BG{a) = Nbg 



u(kFa) 



1 



-= e lka . (15) 



C. Overlap Calculation 

The non-relativistic transition matrix element for a 
transition AA — > H inside a nucleus is given by (sup- 
pressing spin and flavor) 

T { aa } ^h = 2niS(E) J d 3 a d 3 R CM JJ d 3 p i d 3 \ i 

i—a,b 

x M^^e^-^^ (16) 

where 5(E) =5{E H - E AA ), = tfft V ,6 > A a <"), and 
ipnuc = VWc(a) is the relative wavcfunction function of 
the two A's in the nucleus. The notation {AA} is a re- 
minder that the A's are in a nucleus. The plane waves 
of the external particles contain normalization factors 
l/VV and these volume elements cancel with volume 
factors associated with the final and initial phase space 
when calculating decay rates. The integration over the 
center of mass position of the system gives a 3 dimen- 
sional momentum delta function and we can rewrite the 
transition matrix element as 

T {aa} ^h = (2ir) 4 i6 4 (k f - h) M {AA} ^ Hl (17) 

where |A^|{aa}^h is the integral over the remaining in- 
ternal coordinates in (16). In the case of pion or lepton 
emission, plane waves of the emitted particles should be 
included in the integrand. For brevity we use here the 
zero momentum transfer, k — approximation, which 
we have checked holds with good accuracy; this is not 
surprising since typical momenta are < 0.3 GeV. 



Inserting the IK and BBG wavefunctions and perform- 
ing the Gaussian integrals analytically, the overlap of the 
space wave functions becomes 



|Af| A A- 



■H 



2/ 



1 + /2 

x Nbg I 

J T 



3/4 



. 3/2 



d 3 a 



u(k,Fa) 
kpa 



where the factor 1 /4 comes from the probability that two 
nucleons are in a relative s-wave, and / is the previously- 
introduced ratio of nucleon to H radius; an — f olb- 
Since Nbg has dimensions V -1 / 2 the spatial overlap 
A4{aa}— h is a dimensionless quantity, characterized by 
the ratio /, the Isgur-Karl oscillator parameter as, 
and the value of the hard core radius. Fig. 1 shows 
\-^\\ AA }^h versu s the hard-core radius, for a range of 
values of /, using the standard value of as — 0.406 fm -1 



for the IK model[20] and also a B 
parison. 
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Figure 1: Login of \M\ 2 AA ^ H versus hard core radius in 
fm, for ratio / = Rn/Rh and two values of the Isgur-Karl 
oscillator potential, with ag = 0.406 fm -1 (thick lines), and 
a B = 0.221 fm -1 (thin lines). 



IV. WEAK INTERACTION MATRIX 
ELEMENTS 

Transition of a two nucleon system to off-shell AA 
requires two strangeness changing weak reactions. 
Possible AS = 1 sub-processes to consider are a weak 
transition with emission of a pion or lepton pair and an 
internal weak transition. These are illustrated in Fig. 3 
for a three quark system. We estimate the amplitude 
for each of the sub-processes and calculate the overall 
matrix element for transition to the AA system as a 
product of the sub-process amplitudes. 
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Figure 3: Some relevant weak transitions for NN — > HX 



In this way we have roughly estimated all the matrix 
elements for the relevant sub-processes based on weak- 
interaction phenomenology. 



NUCLEAR DECAY RATES 



A. Lifetime of doubly-strange nuclei 



The matrix element for weak pion emission is esti- 
mated from the A — > Nn rate: 



\M\U Nv 



1 2to a 1 



(27r) 4 $ 2 T A ^ NlT 



= 0.8 x 1(T 12 GeV 2 



(19) 

By crossing symmetry this is equal to the de- 
sired \-M\%^ A% , in the approximation of momentum- 
independence which should be valid for the small mo- 
menta in this application. Analogously for lepton pair 
emission we have 



\M\l_ 



1 



2mA 



1 



■Nev 



(2tt) 4 $ 3 T\ — >N CV 



= 3 x 10 



-12 



(20) 



The matrix element for internal conversion, (uds) — > 
(udd), is proportional to the spatial nucleon wave func- 
tion when two quarks are at the same point: 



\M\a-*n ~< ipA\S 3 (ri-r2)\ipN > 



Gf sin 9 C cos 6 C 



m„ 



(21) 



where m q is the quark mass introduced in order to make 
the 4 point vertex amplitude dimcnsionlcss[23]. The ex- 
pectation value of the delta function can be calculated in 
the harmonic oscillator model to be 



< ipA\S 3 (ri-r 2 )\ipN > 



a B 



= 0.4xl0~ 2 GeV 3 



The delta function term can also be inferred phenomeno- 
logically in the following way, as suggested in [23]. The 
Fermi spin-spin interaction has a contact character de- 
pending on o^<72 /to 2 , (5(r*i — r 2 ), and therefore the delta 
function matrix element can be determined in terms of 
electromagnetic or strong hyperfme splitting: 

(m s o - to s+ ) - (m n - to p ) = a^—^ < 6 3 (ri - r 2 ) >(23) 

to a - m N = < S 3 (fi - r 2 ) >(24) 

where m q is the quark mass, taken to be toat/3. Using 
the first form to avoid the issue of scale dependence of 
as leads to a value three times larger than predicted by 
the method used in (22), namely: 

< ipA\S 3 (n - r 2 )\ip N > =1.2xl0~ 2 GeV 3 . (25) 

We average the expectation values (22) and (25) and 
adopt 



\M\l^ N 



4.4 x 10~ 15 . 



(26) 



The decay rate of a doubly-strange nucleus is: 



2(2to aa ) 



(27) 



x <& 2 |-A4| 



AA^H- 



where $2 is the two body final phase space factor, defined 
as in [13], and toaa is the invariant mass of the A's, 
~ 2toa. The factor K contains the transition element 
in spin flavor space. It can be estimated by counting 
the total number of flavor-spin states a uuddss system 
can occupy, and taking K 2 to be the fraction of those 
states which have the correct quantum numbers to form 
the H. That gives K 2 <~ 1/1440, and therefore we write 
K 2 = (1440 K1440)" 1 - Combining these factors we obtain 
the lifetime estimate 



ta 



3(7) Ki44Q lg 
A' it ~ . . 77o 10 S, 



\M\ 2 AA ^ H 



(28) 



where the phase space factor was evaluated for to# = 
1.8(2) GeV. 



(22) „ 




Figure 2: Logio of \M\ AA ^ H versus hard core radius in 



fm, for f=2, 3 and 4 and two values of the Isgur-Karl oscillator 
potential. Thick lines refer to ag = 0.406 fm — 1 , and thin lines to 
a B = 0.221 fin -1 . 



Fig. 2 shows |-M| 2 AA }^ ff in the range of / and hard- 
core radius where its value is in the neighborhood of 
the experimental limits, for the standard choice as — 
0.406 fm^ 1 and comparison value as = 0.221 fm -1 . In 
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order to suppress T(A A a — ► A' H X) sufficiently that some 
A's in a double-A hypernucleus will decay prior to forma- 
tion of an H, we require \M\ AA _ H < 10~ 8 . This is satis- 
fied even for relatively large H, e.g., r H < tn/2.3 (rjv/2.1) 
for a hard-core radius 0.4 (0.5) fm. Thus the observa- 
tion of single A decay products from double-A hyper- 
nuclei cannot be taken to exclude the existence of an H 
with mass below 2toa unless it can be demonstrated that 
r H > 1/2 r N . 



B. Nuclear conversion to an H 

If the H is actually stable (m# < 2m p + 2m e ) two nu- 
cleons in a nucleus may convert to an H and cause nuclei 
to disintegrate. NN — ► HX requires two weak reactions. 
Thus the rate for the process Ann — » A' H irir, i s approx- 
imately 



+A'„7T7r 



K 



(2tt)< 



2(2mAr) 
( \M\ N ^\M\ 



(29) 



AA^H 



V (2m A - m H ) 2 



where the denominator is introduced to correct the di- 
mensions in a way suggested by the AA pole approxi- 
mation. Since other dimensional parameters relevant to 
this process, e.g., m q — mAr/3 or Aqcd, are comparable 
to 2mA — nig and we are only aiming for an order-of- 
magnitude estimate, any of them could equally well be 
used. The lifetime for nuclear disintegration with two 
pion emission is thus 



40 k 



1440 



\M\\ A ^ H 



yr, 



(30) 



taking m# = 1.5 GeV in the phase space factor. For the 
process with one pion emission and an internal conver- 
sion, our rate estimate is 



K 1 



(2-k) 4 



2(2m N ) 

X (\M\N^A7r\M\ N ^A\M\ A A^ H ) 2 

leading to a lifetime, for m# = 1.5 GeV, of 

3 K1440 



(31) 



t Ann—*A' h tt 



\M\l A ^ H 



yr- 



(32) 



If m H ^ 1740 McV, pion emission is kinematically for- 
bidden and the relevant final states are e + v or 7; we now 
calculate these rates. For the transition Ann — » A' H ev, 
the rate is 



K 



(27T) 4 



.a. 



2(2TOAr) 

X (|Al|Ar^A e ^|Al|Ar^A|Al|AA^H) 2 - 



(33) 



In this case, the nuclear lifetime is 

K1440 



T A N 



\M\ 2 



10 5 yr, 



(34) 



AA^H 



taking run = 1.8 GeV. For Ann A' H -f, the rate is 
approximately 



Hl V ' 2{2m N ) 

x $ 2 (|M|^aI^|aa^) 2 , 



(35) 



leading to the lifetime estimate 



T A„ 



2 m 



\M\ 2 



10 b yr, 



(36) 



AA—>H 



for m H = 1.8 GeV. 

One sees from Figure 1 that a lifetime bound of 
> few 10 29 yr is not a very stringent constraint on this 
scenario if to# is large enough that pion final states are 
not allowed. E.g., with K1440 = 1 the rhs of cqn (34) 
is >fcw 10 29 yr, for standard as, a hard core radius 
of 0.45 fm, and tr ~ 1/5 rN - in the middle of the 
range expected based on the glueball analogy. If thh 
is light enough to permit pion production, experimental 
constraints are much more powerful. We therefore con- 
clude that ran ^ 1740 MeV is disfavored and is likely to 
be excluded, depending on how strong limits SuperK can 
give. 



TABLE I: The final particles and momenta for nucleon- 
nucleon transitions to H in nuclei. For the 3-body final states 
marked with *, the momentum given is for the configuration 
with H produced at rest. 



mass 


final state 


final momenta 


partial lifetime 


m H [GeV] 


A' H + 


p [MeV] 


y*K 2 \M\\ K ^ H [yr] 


1.5 


7T 


318 


2 icr 3 


1.5 


-KIT 


170* 


0.03 


1.8 


ev 


48* 


70 


1.8 


1 


96 


2 10 3 



VI. DECAYS OF A QUASI-STABLE H 

If 2m at <rnff < rriN + m A , the H is not stable but it 
proves to be very long lived if its wavefunction is compact 
enough to satisfy the constraints from doubly-strange hy- 
pernuclei discussed in sections II and V A. The limits on 
nuclear stability discussed in the previous section do not 
apply here because nuclear disintegration to an H is not 
kinematically allowed. 
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A. Wavefunction Overlap 



Decay rates and lifetimes 



To calculate the decay rate of the H we start from 
the transition matrix element (16). In contrast to the 
calculation of nuclear conversion rates, the outgoing nu- 
cleons are asymptotically plane waves. Nonetheless, at 
short distances their repulsive interaction suppresses the 
relative wavefunction at short distances much as in a nu- 
cleus. It is instructive to compute the transition ampli- 
tude using two different approximations. First, we treat 
the nucleons as plane waves so the spatial amplitude is: 

Ih^aa = 2m5{E AA -E H ) J \\ d 3 p l d 3 X l d 3 a d 3 R CM 

i—a,b 

x ip H ipT ipl b e i{% » +% "»- %H)RcM . (37) 



The integration over Rcm gives the usual 4D 5 function. 
After performing the remaining integrations leading to 
|-M|jj_>aAj as in (17), the amplitude is: 



\M\ H - 



■AA = 



V 



1 + / 2 



6 /o\ 3/4 / x 3/2 

3 \ I an 



(38) 



o 3 2 2- k N~ k N - 

d 3 a e ~3 Q « a -'^ <1 



The amplitude depends on the size of the H through the 
factor / = tm/th- Note that the normalization Nbg m 
the analogous result (18) which comes from the Bcthc- 
Goldstone wavefunction of A's in a nucleus has been re- 
placed in this calculation by the plane wave normaliza- 
tion factor 1/y/V which cancels with the volume factors 
in the phase space when calculating transition rates. 

Transition rates calculated using eq. (38) provide an 
upper limit on the true rates, because the calculation 
neglects the repulsion of two nucleons at small distances. 
To estimate the effect of the repulsion between nucleons 
we again use the Bethe-Goldstone solution with the hard 
core potential. It has the desired properties of vanishing 
inside the hard core radius and rapidly approaching the 
plane wave solution away from the hard core. As noted 
in section IIIB, Nbg — * l/VV, for a — > oo. Therefore, 
we can write the transition amplitude as in (18), with 
the normalization factor 1/y/V canceled with the phase 
space volume element: 



\M\ H - 



>AA 



2/ 



1 + / 2 



3/4 



\ 3/2 



kFa 



This should give the most realistic estimate of decay 
rates. Table 2 shows the overlap values for a variety of 
choices of rn, hard-core radii, and olb- 



Starting from \M\h->aa we can calculate the rates for 
H decay in various channels, as we did for nuclear conver- 
sion in the previous section. The rate of H — > nn decay 
is 



K 2 



■ (2n) 4 m 5 q 
2 ran 

x (\M\%^ A \M\ H ^ AA y 



<5> 2 {m H ) 



(40) 



The lifetime for this transition, for rnn = 1.9 (2) GeV, is 
tr^nn ~ 5(2) 10 12 K1440 Mo yr, (41) 

where we have introduced \M\ H ^ AA = 710~ 8 /jU ; values 
of Mo > 1 are consistent with hypernuclear constraints. 
Thus we see that the H is stable on cosmological time 
scales if its mass is < 2.04 GeV. 

If 2.04 GeV < m H < 2.23 GeV, H decay requires only 
a single weak interaction, so the rate in eq. (40) must be 
divided by |A^|^-^ A given in eqn (26). Thus we have 



th^na ~ 10 sec K1440 Ho- 



rn 



Finally, if mn > 2.23 GeV, there is no weak interaction 
suppression and 



Tff^AA ~ 10 «i44o Mo sec. 



(43) 



Equations (41)-(43) with Mo = 1 give the lower bound on 
the H lifetime, depending on its mass. 

Our results for the H lifetime are dramatically differ- 
ent from the classic calculation of Donoghue, Golowich, 
and Holstein [24], because we rely on experiment to put 
an upper limit on the wavefunction overlap \M\ 2 h^,aa- 
The bag model is not a particularly good description 
of sizes of hadrons, and in the treatment of [24] the H 
size appears to be fixed implicitly to some value which 
may not be physically realistic. Furthermore, it is hard 
to tell whether the bag model analysis gives a good ac- 
counting of the known hard core repulsion between nu- 
cleons. As our calculation of previous sections shows, 
these are crucial parameters in determining the overlap. 
Our treatment of the color-flavor-spin and weak inter- 
action parts of the matrix elements is rough but should 



TABLE II: |X|i_ AA in GeV~ 3/2 for different values of 1/f 
(rows) and nuclear hard core (columns), for cxb\ = 0.406 fm _1 
and qb2 = 0.221 fm" 1 . 





0.4 fm 


0.5 fm 


OLBl 




OLBl 


OLB2 


0.3 


2 10- 10 


4 10~ 6 


4 10" 13 


6 10" 7 


0.4 


9 10~ 7 


8 10~ 4 


2 10" 8 


3 10" 4 


0.5 


4 10~ 4 


0.02 


1 10" 6 


0.01 
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give the correct order-of-magnitude, so the difference in 
lifetime predictions of the two models indicates that the 
spatial overlap in the bag model is far larger than in our 
model using a standard hard core and taking rjj w 0.4 
fm (which accounts for the suppression measured in hy- 
pernuclear experiments). The calculation of the weak 
interaction and color-flavor-spin matrix elements in ref. 
[24] is more detailed than ours and should be more ac- 
curate. It could be combined with a phenomenological 
approach to the spatial wavcfunction overlap to provide 
a more accurate yet more general analysis. We note that 
due to the small size of the H, the p-wave contribution 
should be negligible. 

One would like to use experiment to place limits on 
the product of the local density of H's and the H decay 
rate, if the H is long-lived enough to be dark matter, 
i.e., ma < ttin + m\. Estimates of the local number 
density of H's at various depths in the Earth, assuming 
the dark matter consists of H's and H's, will be discussed 
in ref. [25]. The Sudbury Neutrino Observatory (SNO) 
can probably place good limits on the rate of H — > nn 
in that detector. The next most important channel H — ► 
77777 should be easy to detect in SuperK for photon energy 
in the low-background range w 20 — 100 MeV [26], or in 
Kamland for lower photon energies. 2 The rate is: 

(2tt)W 

w K chem—z <&z{m H ) (44) 

x (I-MI^aIMI/^aa) 2 

leading to 

t h ^n Ni « 2 10 19 (3 10 17 ) K1440 Mo yr, (45) 

for m H = 1.9 (2) GeV. The lifetime for H — ► npev is 
similar in magnitude. It is more sensitive to ran due to 
the 4-body phase space. 

Z m H 

x (IMIjv^aIMI N^Acv\-M\h^Aa) 

and 

TH^pne. ~ 5 10 19 (3 10 16 ) k 1440 Mo yr. (47) 
for m H = 1.9 (2) GeV . 

VII. SUMMARY 

We have considered the constraints placed on the H 
di-baryon by the stability of nuclei and hypernuclei with 
respect to conversion to an H dibaryon, and we have cal- 
culated the lifetime of the H if it is not stable. We used 



the Isgur-Karl wavefunctions for quarks in baryons and 
the H, and the Bcthc-Goldstone wavcfunction for nuclc- 
ons in a nucleus, to obtain a rough estimate of the H- 
baryon-baryon wavefunction overlap. Observation of A 
decays in double-A hypernuclei is shown not to exclude 
an H as long as < 1/2 rjv- 

Combining our wavefunction overlap estimates with 
phenomenological weak interaction matrix elements, per- 
mits the lifetime of the H and the rate for conversion of 
nuclei to H to be estimated. The results depend radi- 
cally on which channels are kinematically allowed, and 
hence on the H mass, since for each weak interaction 
required there is a substantial suppression in the ma- 
trix element. Our estimates have uncertainties of greater 
than an order of magnitude: the weak interaction ma- 
trix elements are uncertain to a factor of a few, fac- 
tors of order 1 were ignored, a crude statistical estimate 
for the flavor-spin overlap was used, mass scales were 
set to 777at/3, and, most importantly, the models used 
to calculate the wavefunction overlap surely oversimplify 
the physics. The wavefunction overlap is highly uncer- 
tain because it depends on nuclear wavefunctions and 
hadronic dynamics which are not adequately understood 
at present. Nonetheless, the enormous suppression of H 
production and decay rates found in the model calcula- 
tion means that the observation of double-A hypernuclei 
does not exclude ran < 2m a- It is even conceivable that 
an absolutely stable H is possible [5]. 

We calculate the lifetime of the H for various mass 
ranges, taking the H wavefunction to be compact enough 
that hypernuclear constraints are satisfied. If the H 
decays through strong interactions, m# > 2mA, its 
lifetime is > 10~ 9 sec. If its mass is in the range 
777tv + m\ < tuh < 2mA, its lifetime is longer than a few 
10 5 sec, and if 2mjv < m# < mjy + mA the H lifetime is at 
least an order of magnitude greater than the lifetime of 
the Universe. Note that these lifetime bounds do not suf- 
fer from the large uncertainties associated with estimat- 
ing the wavefunction overlap because, for a given value 
of mii, the H lifetime can be related to the hypernuclear 
conversion rate for which we have a rough experimental 
upper limit. It may be that the production rate of H's 
in double A hypernuclei is about equal to the A decay 
rate, so that an appreciable fraction of double A hyper- 
nuclei produce an H. One might hope that in this case 
the H could be observed by reconstructing it through its 
decay products, e.g., H — > Y,~p. Unfortunately, however, 
the long lifetimes implied by the limit on the wavefunc- 
tion overlap mean that the H's would diffuse out of the 
apparatus before decaying. 3 

SuperK can place important constraints on the conjec- 
ture of an absolutely stable H, or conceivably discover ev- 
idence of its existence, through observation of the pion(s), 
positron, or photon produced when two nucleons in an 



2 GRF thanks T. Kajita for discussions on these issues. 



3 GRF thanks K. Imai for enlightening discussions of this topic. 
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oxygen nucleus convert to an H. We estimate that Su- 
perK could achieve a lifetime limit r > few 10 29 yr which 
is in the estimated lifetime range for ran > 1740 MeV 
and ru ~ 1/5 rjy. SuperK and SNO can also place lim- 
its on signatures of H decays if it is not absolutely stable 
yet contributes to the dark matter of the Universe. This 
calculation will be reported elsewhere. The possibility 
that H and anti-H were produced in sufficient abundance 
in the early universe to account for the dark matter and 



baryon asymmetry will also be elaborated elsewhere. 
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